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A 7r-electron network in nanographite systems with zigzag edges exhibits strongly localized edge 
states, which are expected to have peculiar properties. We study effects of electron-electron inter- 
actions on ground-state properties of the systems by means of the weak-coupling renormalization 
group and the density-matrix renormalization-group method. It is shown that the ground state is 
a spin-singlet Mott insulator with finite charge and spin gaps. We also find that the edge states 
are robust against the electronic correlations, resulting in edge effective spins which can flip almost 
freely. The schematic picture for the low-energy physics of the systems is discussed. 

PACS numbers: 71.10.Hf, 73.22.-f, 75.75.+a 



I. INTRODUCTION 

Since the discoveries of fullerenesi and carbon 
nanotubes (CNTs), 2 the nature of "nanographites" , 
graphite-based materials with nano-meter sizes, has at- 
tracted much attention from both fundamental science 
and applications. One of the most striking features of 
the materials is a wide variety of properties they show de- 
pending on their geometrical structure, i.e., size, shape, 
surface condition and so on. For example, CNTs can be 
either metallic or semiconducting depending on the wrap- 
ping superlattice vector^ Such a diversity of possible 
functions realized in nanographites makes them promis- 
ing candidates for nano-scale devices. 

Recently, it has been pointed out that a 7r-electron sys- 
tem in nanographites with zigzag-shaped edges exhibits 
peculiar electronic states strongly localized around the 
edges, which are termed "edge states" i 4 i 5 i 6 Such a state 
does not appear around armchair-shaped edges. The na- 
ture of the edge states has been studied in the graphene 
sheet of nano-meter width, named the "nanographite rib- 
bon" (NGR). Applying a tight-binding model for the 
NCR with zigzag edges, Fujita et al. have shown that the 
system has the electronic states localized on the edges, 
which penetrate from the edges into the bulk decaying 
exponentially—^ The edge states are believed to be re- 
sponsible for a paramagnetic behavior observed in acti- 
vated carbon fibers, which are taken to be an assembly 
of nanographite particles— £ 

Since the edge states are characterized by almost flat 
dispersions in a certain range of the momentum space, 
which result in a sharp peak in the density of state at 
the Fermi energy, possible instabilities to various pertur- 
bations are important and have been investigated to date. 
Concerning the Peierls instability due to electron-phonon 
interactions, it has been shown that no bond alterna- 
tion occurs in the zigzag NGRs with realistic interaction 
strength^ On the other hand, effects of electron-electron 
interactions on the edge states have been investigated by 



a mean-field approximation— iifi The studies have shown 
that an infinitesimal on-site interaction causes a spon- 
taneous spin polarization at the zigzag edges. Similar 
results have been obtained for the zigzag CNTs by the 
density-functional theory— However, it is also known 
that the one-body approximation adopted in these ap- 
proaches is not appropriate for one-dimensional (ID) 
quantum systems such as NGRs and CNTs. In fact, 
Licb's theorem 13 prohibits the spontaneous spin polariza- 
tion for the Hubbard model in the nanographite systems 
with zigzag edges. More detailed analyses with controlled 
approximations are therefore desirable for clarifying the 
low-energy properties of the nanographites. 

In this paper, we study effects of electron-electron in- 
teractions on the nanographite systems with zigzag edges 
using two powerful techniques: the weak-coupling renor- 
malization group (RG) and the density-matrix renormal- 
ization group (DMRG) method. 1 — These approaches al- 
low us to treat the strong quantum fluctuations in a con- 
trolled way for the ID systems such as NGRs and CNTs. 
With the inclusion of electronic correlations, we show 
that the ground state of the zigzag NGR is a gapped 
spin-singlet and the artificially broken symmetry of spin 
rotations in the mean-field type analysis is restored by 
quantum fluctuations. It is rather interesting that the 
edge states survive in the presence of electron-electron 
interactions and play a crucial role in the low-energy 
regime. 

The rest of this paper is organized as follows. The 
model is presented in the next section. In Sec. Ill, we 
show results of the weak-coupling analysis. We first dis- 
cuss properties of the tight-binding model of the NGR 
system. Analytical eigen-wavefunctions and the corre- 
sponding energy spectrum are shown. And then, by 
keeping the low-lying edge states only and taking the 
electron-electron couplings into account, we derive the 
low-energy effective Hamiltonian in the continuum limit. 
It is shown that all coupling terms included in the Hamil- 
tonian are relevant because of a large dynamical exponent 
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FIG. 1: Honeycomb lattice of carbon atoms in the zigzag 
NGRs and CNTs. The open circles and gray squares represent 
carbon sites of sublattice A and B, respectively. 



system size in the x- and y-directions are denoted by L x 
and L y , respectively. Within tight- binding approxima- 
tion, the hopping amplitude t(r,r') is assumed to be t 
between nearest-neighbor sites and otherwise 0. In the 
x-direction, the open boundary condition with L x — ► oo 
is imposed for the zigzag NGR while the periodic bound- 
ary condition is imposed for the zigzag CNT. In the y- 
direction, the open boundary condition is imposed on the 
zigzag edges. We set L y to be even so that the system al- 
ways has the reflection symmetry in the y-direction. All 
carbon sites at the edges are assumed to be terminated 
by hydrogen atoms. 

We note that the system is bipartite and the number 
of sites in each sublattice is equal, N A — Nb- At the 
natural filling, with one 7r-electron per site on average, 
the system is particle-hole symmetric. Thus, one can 
apply Lieb's theorem^ to the system, which prohibits 
spontaneous spin polarizations in the ground state. 



of the dispersion of the edge states. The enhanced den- 
sity of states at the Fermi energy generates finite charge 
and spin gaps, and consequently, the system exhibits a 
spin-singlet ground state without spontaneous symmetry 
breaking. The DMRG results arc presented in Sec. IV. 
We find that in the presence of the Hubbard interaction 
the ground state is a spin singlet with finite charge and 
spin gaps, and no spontaneous spin polarization appears. 
Besides, electrons around the zigzag edges correlate ferro- 
magnetically to compose effective spins which are easily 
polarized by applying an external field, while electrons 
in the bulk form a spin-singlet state which are hardly 
magnetized. The effective spins are localized around the 
zigzag edges even for a rather strong coupling, showing 
that the localization character of the edge states persists 
robustly against the electron-electron interactions. The 
schematic picture for the low-energy physics of the sys- 
tems is discussed. Finally, our results are summarized in 
Sec. V. 



II. MODEL HAMILTONIAN 

In this paper, we consider a 7r-electron system at half- 
filling on a honeycomb lattice with zigzag edges as shown 
in Fig. n The effect of electron-electron couplings is in- 
corporated by introducing the on-site Hubbard interac- 
tion U . The Hamiltonian is 

H = Ho + Hu, (1) 
Ho = -Y, E t(r,r')[ci(r)c a (r')+h.c], (2) 

r,r' a=\,{ 

Hu = Uj2^(r)H(r), (3) 

r 

where c a (r) {c' a (r)) is the fermion annihilation (creation) 
operator at the site r = (x,y) and n a (r) = c< a (r)ca(r). 
The definition of the site index is shown in Fig. ^ The 



III. RENORMALIZATION GROUP 

A. Band structure 

To study the system in the weak-coupling limit U <C t, 
it is natural to consider the case of U = and analyze 
the band structure first. We thus begin with diagonaliz- 
ing the hopping Hamiltonian Ho- Since the Hamiltonian 
is translational invariant along the x-direction, we can 
perform the partial Fourier transformation^ 

c a (k x ;y) = —== ^ exp[-ik x (x + S)]c a (x, y) (4) 

v x X 

where 5 = 1/2 for y = 0, 1 (mod 4) and 6 = for y = 2, 3 
(mod 4). By the transformation, the hopping Hamilto- 
nian Ho is mapped into bond-alternating chains decou- 
pled from each other, 

Ho = Y.Ko(k x ), 

Ly-1 

Ho{k x ) = +1 ) Ca (**5f) + h * c -]' 

(5) 

where t(y) = t\ = 2tcos(k x /2) if y is odd while t(y) = 
t-2 = t if y is even. It is important to note that the 
effective hopping t\ depends on the momentum k Xl and 
consequently, the motions in the x- and y-directions are 
entangled. 

Each of the bond-alternating chains is composed of two 
sublattices A (odd y) and B (even y). The hopping be- 
tween these two sublattices lead to the coupled Harper 
equations, 

- t^(y A ) - t 2 <t>{y A + 2) = Ecj>(y A + l), (6) 
-t 2 4>{yA-l)-ti<ff(y A + l) = E4>{y A ), (7) 
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where y A = 1,3, 5,..., L y — 1 are the lattice points of 
sublattice A. The open boundary condition in the y- 
direction requires (j)(L y + 1) = and (f)(0) = 0. From 
the coupled Harper equations, the wavefunctions of the 
eigenstates and the corresponding energy spectrum can 
be obtained in an analytical form. For the chains with 
t\ > <2 (| Me | < 27r/3), the wavefunction is given by 



<P Py (yA) 
4> Pv {va + 1) 



±sm[p y y A + <p(p y )} 
sm[p y (y A + 1)] 



(8) 

where the extra phase ip, chosen to be in the range < 
fiPy) < 7r /2, depends on the magnitude of momentum 

Pv> 



<p(p y ) = tan 



h-t 2 
h + h 



■ ta,n p y 



(9) 



Finitcncss of L y causes the quantization of p y which sat- 
isfies the following constraint 



(Ly + l)Py + (f{Py) = TO7T, 



(10) 



where m is an integer. The energy spectrum for these 
states is given by 



E(k x ,p y ) = ±Jt\ +tj + 2t 1 t 2 cos(2p y ) 



(11) 



Here we emphasize that the momentum in the y- 
direction, k y , is no longer a good quantum number due to 
the open boundaries. However, the magnitude of the mo- 
mentum p y remains a good one and would be quite help- 
ful for identifying dominating interactions later. Since 
the wavefunction Eq. JSJ is extended into the bulk sites, 
we call the states "extended states" . 

For the chains with t% < *2 i\k x \ > 2tt/3), peculiar 
localized states show up. One finds that there are only 
L y — 2 extended states in Eq. (JHJ). The missing two 
states are linear combinations of localized states near the 
edges with complex momentum k y — ±7r/2 ± ij. The 
wavefunctions for these edge states are 



^(k x ; da) = e 



±sinh[7(L y + 1 - y A )} 
sinh[7(y j4 + 1)] 



(12) 



The imaginary part of the momentum 7 7^ is the inverse 
of localization length of the edge states and satisfies the 
constraint 



(L y + 1)7 = tanh 



ti+h 



tanh 7 



(13) 



It is clear that the nonzero solution of 7 only exists when 
t\ < t%. The energy spectrum for these edge states is 
given by 



E(k x , 1 ) =±^t\ + t 2 2 - 2*1*2 cosh(2 7 ). (14) 

We thereby obtain the whole band structure of the 
tight-binding Hamiltonian, Eqs. and (|14f> . The band 



E(k x ) 




FIG. 2: Band structure of the zigzag NGR with L y = 12. 
The bold lines represent dispersions around the Fermi point, 
k x — TV. 



structure for L y = 12 is shown in Fig. [21 as an exam- 
ple. It is worth noticing that the energy spectrum of the 
edge states is always lower than those of the extended 
states. It means that the edge states would dominate 
the low-energy physics in the weak-coupling limit while 
the extended states in the bulk do not play a crucial role 
because all of them are suppressed by a finite gap of order 

t/Ly. 



B. Renormalization group analysis 

Now, we are ready to treat the interacting Hamilto- 
nian Tin- After the partial Fourier transformation, the 
interacting Hamiltonian is written as, 



U 



^ ' ^ ' fi(k X i k'x 2 k Xi ) 



{*.J|/=1 

x c] (k Xl ; y)c T (k X2 ; y)c\ (k X3 ; y)c l (k Xi ; y). (15) 

The delta function represents the momentum conserva- 
tion in the x-direction. The Hamiltonian can be rewrit- 
ten in terms of the eigenstates of 7io by expanding the 
electron operator on the diagonal basis, 



c a (k x ;y) =y^ y (/>p{y)ip pa (k x ). 



(16) 



where ^ is taken for all the extended states p — p y and, 
if ever, the edge states p = ±7. The resultant Hamilto- 



nian is 



ftu — ^ t ^ t $(k Xl k X2 + k X3 k X4 ) 

X {k Xi }{ Pi } 

y 

X V'p lT (^l)^,, 2 T(^2)^p 3i (fe3)V'p 4 l(^4)-(l7) 
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In general, pi are incommensurate and the dominant con- 
tribution comes from the pairwise equal pi vertices while 
other kinds of vertices only occur in a limited tiny phase 
space. This great reduction of dominant interactions en- 
ables us to classify them into three categories, (i) The 
interactions involving only the edge states: Since the nor- 
malized wavefunction of the edge states is order one near 
the edge and vanishingly small in the bulk, the wave- 
function product of 4> Pi (y) in Eq. i|17|) after summation 
over y is order one. Thus the edge-edge interaction re- 
mains 0(U). (ii) The interactions involving a pair of 
the extended states and a pair of the edge states: Due 
to the spreading of wavefunction in the extended state, 
the summation of the wavefunction product in Eq. i(T7|) 
is roughly |^ Ph (1)| 2 ~ 1/L y . The edge-bulk interaction 
is therefore only 0(U/L y ) and becomes smaller as the 
width of the NGR grows, (iii) The interactions involv- 
ing only the extended states: From the similar argument 
to the case of the edge-bulk interaction, the bulk-bulk 
interaction turns out to be 0(U). The effect of this in- 
teraction in the bulk system (L y — > oo) has been studied 
in Refs. and by the weak-coupling analysis. It 
has turned out that the interaction opens a finite energy 
gap even if the noninteracting Hamiltonian H.q includes 
gapless chains with k x — ±27r/3. 

Let us move on to deriving the low-energy effective 
Hamiltonian in the continuum limit. In the weak- 
coupling limit, the effective theory is dominated by the 
low-lying edge states near k x sa 7r because all other bands 
in the bulk are gapped. Thus, we may keep only the edge 
states and neglect all the gapped modes. That is to say, 
only the kinetic energy of the edge modes and the edge- 
edge interaction term must be considered in the limit. It 
is useful to construct the antibonding and bonding fields 
ip± a (x) from the field for the edge states defined in Eq. 
(El), 



A dk 
-A 2tt 

1 



ip± a {k x +n), (18) 



1p±a(k x ) = -j=[lp +ia (k x ) ± V>- 7 a(M]> ( 19 ) 

where A is a momentum cutoff. Furthermore, we approx- 
imate the energy spectrum near k x « 7r to the lowest 
non- vanishing order E±(k) ~ ±c(k x — ir) z , where the ± 
signs stand for the antibonding and bonding edge modes 
and the dynamical exponent z = L y — 1. The kinetic en- 
ergy can then be expressed in terms of the fields ip± a (x) 
in the continuum limit 



Ho = dx 



-± 



Pi>l a (x)(-id x y^ pa ( x ), 



(20) 



where the constant c is set to unity for notation simplic- 
ity. 

While it is rather straightforward to write down the 
kinetic terms in the form of Eq. (|20|) , expressing the edge- 
edge interaction in terms of the field operators is far from 
trivial. Our strategy to go from lattice to continuum limit 



is to express the lattice operator for electrons in the edge 
states, c a (x, y), in terms of the field operator tp± a (k) and 
relate it to 4'±a(x) by Fourier transformation. To achieve 
the goal, we need to express c a (x,y) in the eigenbasis of 

Tio, 



c a (x,y) = J T^e lfe 



p=± 



4* P {y)^pu{k). 



(21) 



Here a difficulty comes from the entanglement of summa- 
tion over k x and P. Since 7 depends on the longitudinal 
momentum k x as in Eq. Ijl3|l . the wavefunction 4> P {y) 
also has an implicit dependence on k x . As a result, the 
lattice operator c a (x, y) is not simply related to the con- 
tinuous field ip±a(k) by Fourier transformation which re- 
quires the change of order between the summation over P 
and the integral over k x . However, since the edge states 
with k x ~ 7r are localized near the edges, it is reasonable 
to adapt the "sharp-edge" approximation which assumes 
that all the wavefunctions of the edge states take the 
simple form as for k x = tt, i.e., 



^±{y) = -/|(<W ±( W<y)- 



(22) 



Within the sharp-edge approximation, the relation be- 
tween the lattice and the field operators becomes very 
simple, 

CaCM) = -^[lp +a {x)+^_ a (x)\, (23) 
e iirx 

~c a (x,L y ) = —[^ +a { x )-^_ a {x)\, (24) 

c a {x,y) = (2 < y < L y - 1). (25) 

By substituting Eqs. f2l^l - lt2*5|) into the edge-edge in- 
teracting Hamiltonian 



7i e 



U ^2c\(x,y)c^(x 1 y)c\(x,y)c l (x,y), (26) 



we obtain its continuum counterpart, 

iiee = 9c{J\ + J-) + 9pJ+J- - 9aJ+ ' J- 

+ |(7 + 7l + /_4). 
Here the SU (2) invariant currents are defined as 

Jp = 2^«^p«' 

Jp = Tj^iaVap^pf}, 
Ip = ^i a e a0tppf3, 



(27) 

(28) 
(29) 
(30) 



where <x and e a p denotes Pauli matrices and an an- 
tisymmetric tensor, respectively. The bare values of 
these couplings generated by the on-site interaction U 
are g c = g p = g a = g u — U. We note that Eq. l(77|) 
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is the most general form of the interacting Hamiltonian 
which includes four-fermion interactions with preserving 
the U(l) x SU(2) symmetry and the momentum conser- 
vation. In fact, one can use Eq. I|27l) to consider more 
complicated short-range interaction, which would gener- 
ate a different set of bare couplings. 

Finally, the low-energy field theory of the zigzag NGR 
in the continuum limit is described by Eqs. i|20|) and (|27ll . 
At the tree level, the scaling dimensions of the all four- 
fermion interactions are d{g) = z — 1 by simple dimen- 
sion counting. Since the dynamical exponent is greater 
than one for the NGR, all the couplings g c ,gp,g<?, and 
g u are relevant under the RG transformation. However, 
we note that relevance of these couplings does not nec- 
essarily warrant excitation gaps in the energy spectrum. 
To identify the ground state, it is necessary to apply an- 
other approach such as the bosonization technique. We 
will discuss it in the next subsection. 



C. Mott insulating ground state 

In this subsection, we consider what we can derive from 
the low-energy effective Hamiltonian obtained in the pre- 
vious subsection. One clue is the fact that the effective 
Hamiltonian has a form strikingly similar to the usual 
effective theory for ID systems except for the large dy- 
namical exponent z > 1. Hence, we discuss the analogy 
between the NGRs and a one-chain system to derive the 
ground-state properties of the NGR. On the other hand, 
the spectrum of the low-lying edge states in the NGR 
can be viewed as a two-chain system, whose bonding and 
antibonding bands touch the Fermi energy at one Fermi 
point as discussed below. Therefore, it is also useful to 
consider the NGR in the analogy with the two-chain sys- 
tem. We will see that the both analyses give the same 
conclusion: the NGR has a spin-singlet ground state with 
finite charge and spin gaps. 

Let us start with the first approach, i.e., considering 
the NGR in the analogy to a one-chain system at half- 
filling. In this analogy, the field ip +a (ip- a ) is regarded as 
a right- (left-) moving field operator. Since the dispersion 
of the fields ip±a{x) is not linear, the rigorous mapping 
from the fermionic theory to the bosonic one is not possi- 
ble. However, it is likely that the large dynamical expo- 
nent z just serves to make all interactions relevant and we 
can still use the bosonization techniques to clarify physi- 
cal properties of the ground states qualitatively. That is 
to say, we can treat the dynamical exponent z = 1+e as a 
perturbation and still use the conventional bosonization 
rules to study the physical properties of the NGRs. Un- 
der the assumption, we follow conventional bosonization 
steps and introduce two pairs of bosonic fields 

ip± a (x) ~ VP e^^ ±e ^/ 2 , (31) 

where <fi a (x) and a (x) are respectively the phase and 
displacement fluctuations, which obey the commutation 



phase 


(V20 CTJ V26 p ) 


C s 


D 


New 


(tt.O) 








CDW 


(O.tt) 







Dimer 


(0,0) 








TABLE I: Comparison between different phases, 
relation, 

[<f> a (x),6 a ,(x')} = -2Tri[l + sign{x-x')]6 aa ,, (32) 

Furthermore, we switch to the charge and spin basis de- 
fined by 

(^ CT ,^ ff ) = -i=[(0 T ,0 T )±(^,0 i )]. (33) 

The effective interaction which is relevant to gap forma- 
tion in Eq. I|27|l can be written down in terms of these 
bosonic fields, 

H ee = g a cos(\/26» CT ) - g u cos(V26» p ). (34) 

The absence of their dual fields 4> p and <^ CT in the inter- 
actions comes from the charge and spin conservations. 
For the on-site interaction, both g a and g u are positive 
and grow under RG transformation. Consequently, the 
corresponding bosonic fields arc pinned at specific values 
up to fluctuations, 

V26 a ~ 7T, V29 p ~ 0. (35) 

The signs of the couplings are of crucial importance here. 
Expanding the cosine term around the pinned values, the 
system acquires finite gaps in both charge and spin sec- 
tors. Therefore, the ground state is a Mott insulator with 
both charge and spin gaps. The formation of the charge 
gap is not surprising because the electron density is com- 
mensurate with the underlying lattice structure. The 
finite spin gap indicates that the spin rotational symme- 
try is not broken spontaneously and the ground state is a 
spin singlet, as predicted by Lieb's theorem. The result 
differs from the mean-field theory predictions where the 
quantum fluctuations are ignored. We believe that the 
strong fluctuations in ID systems would invalidate the 
mean-field approach and restore the symmetry which is 
broken at the mean-field level. 

It is worth noting that for the usual weakly-coupled ID 
system with linear energy spectrum, the bosonic fields 
never get pinned at the values listed in Eq. I|35l) because 
it is not stable in the RG flow. However, the new phase 
here is stable due to the large dynamical exponent which 
corresponds to the van Hove singularity in the density of 
states. We summarize the pinned values for each phase in 
Table [Q To differentiate this new phase from the famil- 
iar charge density wave or the dimer phase, one needs to 
compute the order parameters, i.e., charge density mod- 
ulation C s and bond dimerization D defined as 

C s = (-ir(ai(x)a a (x)) 
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- (^-a) + <^ + a>. (36) 

D = { -^-(ai( X )a a (x + l) + aUx + l)a a (x)} 



(37) 



where a a (x) is the lattice operator in the one-chain sys- 
tem and related to the fields as a a (x) ~ e ikpx ip +a (x) + 
e~ tkFX ip„ a (x) with k F = 7r/2. Using the relations 



ie -ie p /y/2 e -He„/V2 



(38) 
(39) 



and substituting the pinned values of the bosonic fields, 
one obtain the expectation values of the order parame- 
ters. The results are listed in Table [I] It is clear that 
there is no long-range order or symmetry breaking for 
the ground state of the zigzag NGRs. 

Now we take another perspective and view the NGR as 
a system of two chains coupled via an interchain hopping. 
Figure E3 (a) shows the band structure of the two-chain 
system for the special point t Tuns /ti eg — 2, where i rung 
(iieg) is an amplitude of the interchain (intrachain) hop- 
ping. As the hoppings are tuned at the special value, 
the chemical potential cut through the bonding (anti- 
bonding) band at single Fermi point k x — tt {k x = 0) 
with vanishing Fermi velocity v F — > 0, rather than two 
pairs of Fermi points in general. On the other hand, both 
the bonding and antibonding bands of the low-lying edge 
states in the NGRs have single Fermi point k x = tt with 
v F — > 0. Hence, in the low-energy limit, the peculiar 
band structure in the NGRs can be viewed as a two-chain 
system with shifting the momentum of the antibonding 
field k x -> k x + tt [See Fig. H and Fig. (b)]. One 
might worry that the dynamical exponent in the NGRs 
z = L y — 1 may be different from that in the two-chain 
system z — 2. However, the dynamical exponent z in 
this limit only affects how the Fermi velocity approaches 
zero and should not give rise to qualitative difference. 

In the weak-coupling limit, the lattice electron opera- 
tor diotix) in the i-th chain of the two-chain system can 
be described by two pairs of chiral fields 

p i-KX 

a la (x) = -^{^ +a {x)+e^ x ^ a {x% (40) 

a 2a (x) = -j=-[xj} +a {x) - e lQx ij)- a (x% (41) 

where Q = tt is the momentum difference between Fermi 
points of bonding and antibonding bands in the two-chain 
system. The relation of lattice operators to the field op- 
erators for the two-chain system in Eqs. I|40|l and i|41|) 
are strikingly similar to those for the NGR in Eqs. lll'-'il) 
and (|24p. except the momentum shift Q. This enables us 
to translate the excitations and correlations in the two- 
chain system directly into those for the NGRs. 

For on-site repulsive interaction, the ground state of 
the two-chain system is the Mott insulator with finite 




% K 



FIG. 3: Schematic picture of the band structure of (a) the 
two-chain system with t luns /ti cg — > 2, and (b) the same sys- 
tem with the momentum shift Q = n for the antibonding field. 
The bold lines represent dispersions around Fermi points. 



charge and spin gaps. In addition, there exists neither 
broken symmetry nor local order parameter. This Mott- 
insulating phase resembles the peculiar phase we dis- 
cussed above by making analogy to the one-chain case. 
Thus, both approaches support that the ground state of 
the NGRs is a spin singlet with finite charge and spin 
gaps. 

The two-chain analogy also gives us information about 
the spin-spin correlation. It is known that the two-chain 
system has magnon excitations with spin S — 1 and 
momentum P = (tt,tt). The momentum shift Q in the 
x— direction shifts the low- lying magnon to P = (0, tt). 
This means that the spin-spin correlation functions in 
the NGRs are ferromagnetic within the same edge and 
antifcrromagnetic between opposite edges. This result is 
consistent with the fact that the sites in the same edge be- 
long to the same sublattice while those in opposite edges 
belong to different sublattices. 



IV. NUMERICAL RESULTS 

In this section, we present numerical results for the tt- 
electron system Q at half- filling. Hereafter, we set t = 1. 
Using the DMRG method with improved algorithmji2ii£ 
we have calculated various energy gaps and correlation 
functions. The number of kept states is up to m = 1200 
per block. Numerical errors due to the truncation of the 
DMRG calculation are estimated from the difference be- 
tween the data with different m's. We discuss the results 
for the zigzag NGRs and the zigzag CNTs in the following 
subsections. 



A. Zigzag nanographite ribbons 



In this subsection, we show results on the zigzag NGRs. 
We first discuss the results for L y — 4. In all the cases 
calculated, we have found that the ground state belongs 
to the subspace of M = 0, where M is the total magneti- 
zation in the system. In Fig. 21 we show the data of the 
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FIG. 4: Lz-dependence of the energy gaps in the NGRs 
for U — 1; (a) the charge gap A c , (b) the spin gap As . 
The solid and open circles represent the data for L y — 4 and 
L y — 6, respectively. The numerical errors due to the DMRG 
truncation are smaller than symbols. The gaps for L y = 4 
are extrapolated by fitting the data to a polynomial form, 
A(L X ) = A(oo) + a/L x + b/L x where a and b are fitting pa- 
rameters, (c) The [/-dependence of the extrapolated gaps for 
K = 4. 



charge and. spin gaps defined as 




where E§{N^,N\) is the lowest energy in the subspace 
of JV-f up-spin and iVj down-spin electrons and N is the 
total number of sites. For L y — 4, both the energy gaps 
are extrapolated to non-zero values at L x — ► oo. The in- 
dependence of the extrapolated gaps suggests that both 
of the gaps start to open from infinitesimal U , i.e., U c = 
0. We thus conclude that for arbitrary U > the ground 
state of the zigzag NGRs is a spin singlet with finite 
charge and spin gaps, being consistent with the weak- 
coupling analysis. 

To see how the spin polarizations are distributed in 
real space, we have calculated the local spin polarization 

(S z (r)) M = l(ni(r)- H (r))M, (44) 

where (• • •) m denotes the expectation value in the lowest- 
energy state in the subspace of M. We have found that in 
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(a) L y = 4, U = 1 




(b) L y = 6, U = 1 




FIG. 5: Distribution of local spin polarization {S z (r))M=i 
in the NGRs with (a) L y = 4 and (b) L y = 6. The Hubbard 
coupling and the total magnetization are U — 1 and M = 1. 
Open and gray circles respectively represent the positive and 
negative values of the spin polarization while the areas of 
them are proportional to the magnitude. 

the ground state the local spin polarization (S z (r))M=o 
is zero at all sites within the numerical error, supporting 
that the ground state is spin singlet. There is no spon- 
taneous spin polarization. Figure [S] (a) shows the distri- 
bution of the local spin polarization in the magnetized 
state of M = 1 for U = 1 and L y = 4. As seen in the fig- 
ure, the magnetization appears almost only at the edge 
sites. This suggests that the electrons in the edge states 
are polarized by an external field more easily than those 
in the bulk. We note that the localization character of 
the magnetization seems to persist even for a rather large 
coupling [7 = 4. In other words, the edge states survive 
robustly against the electron-electron interactions even if 
U becomes larger than the hopping amplitude t. 

In Fig|Sl we show the spin-spin correlation functions 
in the edge sites, 

(S z (x,l)S z (x',l)) M=0 . (45) 

To reduce boundary effects, the correlations are calcu- 
lated between the sites of the same distance from the 
center of the NGR. The correlations turn out to be al- 
ways positive, suggesting a ferromagnetic coupling be- 
tween the sites in the same edge. We have also found that 
the spin correlation between the sites in opposite edges is 
always negative, i.e., antifcrromagnctic. These behaviors 
are consistent with the expectation from the two-chain 
analogy discussed above. The correlations decay expo- 
nentially for all U's reflecting the finite spin gap^ and 
are enhanced as U increases. 

Finally, we discuss the dependence of low-energy prop- 
erties of the zigzag NGRs on the width L y . Since the 
low-energy physics of the NGRs is dominated by the 
edge states which are strongly localized around the zigzag 
edges, it is likely that the singlet-triplet spin gap A^, 
which corresponds to the energy scale required to mag- 
netize the edge spins, comes from the antiferromagnetic 
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FIG. 6: Ground-state spin-spin correlation functions in the 
NGRs with L y = 4 and L x = 20. The circle, square, and tri- 
angle represent the data for U = 1, 2, and 4, respectively. The 
error bars represent the numerical error due to the DMRG 
truncation. 

effective coupling between the electrons in the edge states 
of opposite sides. Thus, it is natural to expect that the 
spin gap decreases as the width L y increases. To see 
this, we have performed calculations for the NGR with 
L y = 6i2& As can be seen in Fig. (b), the spin polar- 
izations in the magnetized state of M — 1 are found to 
be strongly localized at the zigzag edges, suggesting that 
the edge states exist also for L y = 6. Furthermore, we 
can clearly see in Fig. 0]that the spin gap for the same 
L x decreases as the width of the ribbon L y increases, as 
expected. Unfortunately, the system size calculated for 
L y = 6 is not large enough to obtain the extrapolated 

gap As 1 ' 1 (La; — > oo). Calculations for larger systems are 
desirable to clarify the situation for the NGRs with larger 



B. Zigzag nanotubes with open edges 

For the zigzag CNTs with open ends, we can apply 
basically the same argument as that for the zigzag NGRs. 
The difference between the two systems is merely the fact 
that in the CNTs the momentum along the x-direction is 
quantized as k x — 2nn/ ' L x due to the finite circumference 
L x . To see what happens, we have done calculations 
for the zigzag CNT with L x = 2. In this case, there is 
one edge state at each zigzag open end. Although the 
width of the CNT is quite narrow, we will show that 
characteristic magnetic properties expected for general 
zigzag CNTs already appear. 

We have found that the ground state always belongs to 
the subspace of M — for all values of U > 0. In Fig. [7| 
we show the data of the charge gap A c and the spin gaps 
As 1 '' and As 2 '. It can be seen in the figure that all the 
gaps take nonzero values in all the case of U > and L y 
calculated. We therefore conclude that the ground state 
of the zigzag CNTs with finite L y is a spin-singlet Mott 
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FIG. 7: Charge gap A c and spin gaps As 1 ' and As ' of the 
CNTs with L x = 2 for U = 1,2, and 4 as functions of the 
length of CNTs L y . The numerical errors due to the DMRG 
truncation are smaller than symbols. The dotted lines are 
guides for eye. 



insulator, same as the NGRs. An interesting observation 
here is that the spin gap Ag 1 ' decays exponentially as L y 

increases while A^ 2 ' converges to finite values at L y — > oo. 
Hence, for large enough L y , the magnetization M = 1 can 
be induced by applying an infinitesimal field whereas a 
finite field is needed to magnetize the system to M > 2. 

Next, we discuss how the spin polarizations are dis- 
tributed in real space. Same as the NGRs, the local 
spin polarization in the ground state is zero at all sites. 
Figure [8] shows the distribution of the local spin polar- 
ization (S z (r))M in magnetized states for U = 1. The 
data clearly show that in the state of M = 1 the in- 
duced magnetization are strongly localized at the ends of 
the zigzag CNT. This demonstrates the robustness of the 
localization character of the edge states against the cou- 
pling U. In the state of M = 2, on the other hand, the 
distribution of (S z (r)} m=i seems to be a superposition 
of the edge magnetization and magnetization excited in 
the bulk sites. We have found similar distributions of 
(S z (r)) M for U = 2 and 4. 

The results on the spin gaps and the local spin po- 
larizations lead us to the conclusion that there are two 
different energy scales in magnetic excitations in zigzag 
CNTs: the energy to magnetize electrons in the edge 
states and the one to magnetize the bulk electrons. As 
the length of the CNT increases, the former decays ex- 
ponentially while the latter converges to a finite value, 
the spin gap of the bulk system. This suggests that one 
can regard the edge states and the bulk electrons as two 
almost independent objects; The edge states are essen- 
tially decoupled from the bulk singlet state in spite of the 
presence of U ~ 0(t). Very recently, it has been pointed 
out from a weak-coupling analysis that thicker CNTs in- 
cluding more than one edge states at each open ends also 
exhibit such a decoupling of the edge and bulk statesi^i 
We therefore expect that the conclusion above is valid 
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(a) M = 1 (b) M=2 




circumference circumference 



FIG. 8: Distribution of local spin polarization in the CNTs 
with L x = 2 for U — 1 and L y = 14. The total magnetization 
is (a) M = 1 and (b) M = 2. Open and gray circles respec- 
tively represent the positive and negative values of the spin 
polarization while the areas of them are proportional to the 
magnitude. 

for general CNTs with zigzag open edges. 

C. Picture of the low-energy states 

From the above results, we can deduce a schematic 
picture to represent the low-energy physics of the zigzag 
NGRs and CNTs. In this picture, the 7r-electron system 
consists of two parts: the electrons in the bulk forming 
a spin-singlet state and those in the edge states at each 
zigzag edge. The electrons in the same edge are cor- 
related ferromagnctically with each other to compose a 
large effective spin. The effective spins interact via an ef- 
fective antiferromagnetic coupling across the bulk singlet 
state. Hence, the ground state of the system is a spin sin- 
glet in total. The effective coupling becomes smaller as 
the distance between the effective spins becomes larger, 
and finally, when the distance becomes large enough, the 
effective spins can flip freely, giving a paramagnetic re- 
sponse. The bulk electrons remain forming a spin singlet 
unless a field corresponding to the bulk spin gap is ap- 
plied. We note that the picture above is basically consis- 
tent with the one discussed in Refs. ^3 andl22l 

V. DISCUSSION 

In this paper, we have studied low-energy properties of 
nanographite systems with zigzag edges in the presence of 



on-site Hubbard interactions U using the weak-coupling 
RG and the DMRG method. We have analyzed the Hub- 
bard model on the zigzag NGRs and the zigzag CNTs. 
We find that in both systems the ground state for U > 
is a spin-singlet Mott insulator with finite charge and 
spin gaps. It is also found that the localization property 
of the edge state persists even for a rather large value 
of U, resulting in the effective spins localized around the 
zigzag edges. 

Finally we wish to touch upon further extensions of 
this study. One issue to be studied is effects of fur- 
ther hoppings and long-range electron-electron interac- 
tions. As for magnetism, these long-range terms violate 
the assumptions of Lieb's theorem and open the possi- 
bility of exotic magnetic states. However, it is also ex- 
pected that the spin gap observed in the present study 
tends to stabilize the spin-singlet ground state against 
perturbations. Very recently, it has been pointed out by 
the density-functional theorjii^ that zigzag CNTs with a 
finite length can exhibit a high-spin ground state depend- 
ing on the circumference L x , although the method may 
underestimate effects of strong quantum fluctuation in 
the nanographite systems. To clarify the situation, anal- 
yses with controlled approximations are needed. Another 
issue is stacking effects. Nanographite materials often 
take a stacking form, e.g., nanographite particles com- 
posed by a few graphene sheets and multiwall CNTs. 
Hence, it is important for applications to study effects 
of interlayer interactions. For nanographite particles, it 
has been shown that their magnetic properties change 
drastically depending on types of stacking geometry^ 
In multiwall CNTs, the difference of radius and chiral- 
ity between inner and outer layer is expected to affect 
not only the edge states but also the bulk properties in a 
complicated way. Further extensive studies are desirable 
to tackle the issue. 
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